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Using angle-resolved photoemission (ARPES), it is revealed that the low-energy electronic ex- 
citation spectra of highly underdoped superconducting and non-superconducting La2-xSr cc Cu04 
cuprates are gapped along the entire underlying Fermi surface at low temperatures. We show how 
the gap function evolves to a d x i_ y i form as increasing temperature or doping, consistent with the 
vast majority of ARPES studies of cuprates. Our results provide essential information for uncovering 
the symmetry of the order parameter(s) in strongly underdoped cuprates, which is a prerequisite for 
understanding the pairing mechanism and how superconductivity emerges from a Mott insulator. 
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In the Bardeen-Cooper-Schriefer (BCS) theory of su- 
perconductivity, the symmetry of the superconducting 
gap reflects the order parameter of the superfluid state 
and is directly tied to the symmetry of the interac- 
tions driving the formation of Cooper pairs. Similarly, 
other ordered phases, such as charge- or spin-density 
wave states, can induce gaps whose symmetries are con- 
nected to the underlying order parameters. Thus in 
high-temperature superconductors, where superconduc- 
tivity is found in close proximity to magnetic and charge 
order, the symmetry of the gap function is of critical 
theoretical importance. It is now widely accepted that 
the superconducting gap in moderately hole-doped high- 
temperature superconducting copper oxides (cuprates) 
exhibits a node located along the diagonal (0,0)-(7r,7r) 
line of the Brillouin zone (BZ) [TJ |2], consistent with 
an overall gap function of pure d x i_ y i symmetry [3]. A 
key issue is whether d x 2_ y 2 is the only form of the gap 
function for all cuprates over the full range of dopings. 
Using angle-resolved photoemission (ARPES) we reveal 
that the low-energy electronic excitation spectra of highly 
underdoped superconducting and non-superconducting 
La2- x Sr a: Cu04 (LSCO) are gapped along the entire un- 
derlying Fermi surface (FS) at low temperatures. On the 
zone diagonal, gapless excitations appear as the temper- 
ature and/or doping is increased, and the gap function 
evolves to a d x i_ y i form. 

High quality single crystals of superconducting and 
non-superconducting La2_ a; Sr :c Cu04 were grown using 
the traveling solvent floating zone method. ARPES ex- 
periments were carried out at the Surface and Interface 
Spectroscopy beamline at the Swiss Light Source (SLS). 
Circularly polarized light with hv = 55 eV was used. 
The spectra were recorded with a Scienta R4000 electron 
analyzer. The energy resolutions was about 14-17 meV. 

In Fig. [TJa)-(d) we show ARPES spectra below and 
above the superconducting transition temperature (T c ) 
for highly underdoped superconducting LSCO (x = 0.08, 
T c = 20 K) along the diagonal line of the BZ. The spectra 



were obtained by deconvoluting the raw ARPES data to 
remove the broadening due to the finite instrumental res- 
olution and then dividing the deconvoluted spectra by a 
Fermi distribution function [deconvolution - Fermi func- 
tion division (DFD) method] [4 . Relative to Ep, a gap 
is clearly observed both below T c (10 K) and above T c 
(54 K). The gap closes above ~ 88 K. To reveal the de- 
tails of the gap, we trace the dispersion in the vicinity of 
E F . In Fig. [jje)-(h) we plot energy distribution curves 
(EDCs) from Fig. [TJa)-(d) along the zone diagonal cut. 
At low temperatures (10 K and 54 K), moving from (0,0) 
to (7r,7r), the peak position of the EDCs approaches Ep, 
but before reaching Ep it recedes to higher binding en- 
ergies (Fig. nTe)-(f)). The spectral peak reaches Ep at 
88 K (Fig. fTTg)) and crosses Ep at higher temperatures 
[Fig. [ljh)], signaling the closure of the gap. In Fig. [I] 
(i)-(j) the dispersion bending back below Ep at 10 K and 
crossing Ep at 137 K are better visualized by normalizing 
each EDC to its peak intensity. 

In contrast to the observation of an energy gap along 
(0,0)-(7r,7r) at low temperatures for LSCO (x = 0.08), 
we did not find such a gap in the spectra of the opti- 
mally doped LSCO (x = 0.145, T c = 33 K) sample. The 
spectral peak in this higher hole-doped sample continu- 
ously moves to higher energies and crosses Ep, even at 
temperatures down to 12 K [Fig. [TJl)-(o)]. 

A general observation in cuprates is that when T c 
is crossed by decreasing the temperature, EDC peaks 
near the Fermi momentum (kp) along an off-nodal cut 
sharpen, an energy gap opens at hp, and the renormal- 
ization of the band dispersion becomes more pronounced. 
This behavior was also observed in our ARPES spec- 
tra for LSCO (x = 0.08) along the diagonal cut [red 
arrow in Fig. [I] (k)]. In Fig. [2ja)-(b) we plot the 
EDCs at the kp on the zone diagonal line as a function 
of temperature. The EDC peak width at 10 K (< T c ) 
is considerably smaller than when the spectrum is mea- 
sured above T c (T = 54 K) [Fig. |2jb)]. At low tem- 
perature (10 K) the peak width of the superconducting 
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sample LSCO (x = 0.08) is also smaller than that of 
the non-superconducting sample (x = 0.03) [Fig. [2^b)] , 
which indicates the change in the peak width is associ- 
ated with the superconducting transition. The anoma- 
lous "kink" in the dispersion along the nodal direction, 
which is well-known from previous ARPES studies [5H7] , 
becomes more pronounced below T c [Fig. [5|c)], showing 
that the renormalization is enhanced in the supercon- 
ducting state. To gain more insight about the differences 
between the spectra in the superconducting and non- 
superconducting phases, Fig. [2^d)-(f ) shows the EDCs 
as a function of doping at 10 K. The EDC peak posi- 
tions of the x = 0.03 and 0.08 samples occur at higher 
binding energy than in the x = 0.105 and 0.145 samples 
[Fig. |2^d)], in which a simple d x 2_ y 2 superconducting gap 
was observed [8j [9] . An energy gap of ~ 20 meV along 
the (0,0)-(7r,7r) is also seen in the non-superconducting 
sample (x = 0.03) at 10 K. For the superconducting 
samples, after aligning the peaks to the same position, 
the falling edges of the EDC peaks are almost identical 
at low binding energies [Fig. [2|e)], which demonstrates 
that the coherent peaks have the same width in the su- 
perconducting state, and this width is smaller than that 
of the non-superconducting sample. We observe that un- 
derdoping increases the transfer of spectral weight from 
the coherent peak to high binding energies [Fig. |2^f)]. 

The momentum dependence of the gap as a function 
of temperature is shown in Fig. [3J The measured raw 
EDCs at kp were symmetrized to remove the effects of 
the Fermi function [10]. In the superconducting state, 
the extracted gap sizes were defined as half the peak-to- 
peak separation of the symmetrized EDCs. Above T c , for 
those spectra having no coherent peak, the gap was de- 
fined as half the distance between the two locations where 
the slope has the largest change, as indicated by vertical 
lines in Fig. ^c). For LSCO (x = 0.08), the energy gap 
is highly anisotropic [Fig. [3^f)]. It has a maximal value 



FIG. 1: (Color online) ARPES spectra for 
LSCO with x = 0.08 (T c = 20 K) and 
x = 0.145 (T c = 33 K). (a)-(d) Intensities 
along the zone diagonal [red arrow in (k)] 
at T = 10 K, 54 K, 88 K and 137 K for 
x = 0.08. The spectra were obtained by 
the DFD method 4 . (e)-(h) EDCs from 
(a)-(d) in the vicinity of Jzf- (i)-(j) Images 
in the vicinity of Jsf at 10 K and 137 K. 
Each EDC is normalized to the intensity 
at the peak position (red circles), (k) The 
first quadrant of the BZ. The red curve is 
the FS of LSCO x = 0.08 obtained from 
tight-binding (TB) fits to the experimen- 
tal data. The red arrow indicates the cut 
along which the ARPES data were taken. 
(l)-(m) the same as (a)-(d) but for LSCO 
(x = 0.145) at T = 12 K and 40 K. (n)-(o) 
EDCs from (l)-(m) close to 



at the zone boundary (0 = 0°) and decreases monoton- 
ically along the FS to a minimum at the zone diagonal 
(0 = 45°). However, at low temperatures, the gap (A) 
function strongly deviates from a pure d x 2_ y 2 form, 

A d x2 _ y2 (fe) = A° dx2 _ y2 [cos(k x a) - cos(^a)]/2, (1) 

where a is the lattice constant. Below T c at the hp on 
zone diagonal, a finite gap (hereafter we use the term 
"diagonal gap" for simplicity) is observed, which has an 
amplitude of ~ 20 meV. As the temperature is increased, 
the diagonal gap monotonically decreases. It disappears 
at ~ 88 K, at which point the gap function is very close 
to a pure d x 2_ y 2 form [Fig. |3^f)]. At higher tempera- 
tures, gapless excitations appear on a portion of the FS 
centered at the zone diagonal (Fermi arc) and the arc 
length increases with temperature, a phenomenon that 
has been observed in early ARPES work of underdoped 
Bi2Sr2CaCu208+5(Bi-2212) [11] and other cuprates. In 
contrast to the temperature dependence of the diagonal 
gap, the gap size near the zone boundary is insensitive to 
temperature in the range of 10 K - 150 K; it always has 
a value of - 36 meV [Fig. gf)]. 

The gap function obtained from the symmetrization 
method, in which particle-hole symmetry is assumed, is 
confirmed with an independent procedure. We have col- 
lected high-statistics ARPES spectra [Fig. [IJa)-(f)] along 
some selected cuts [Fig. Qh)] and then applied the DFD 
method to trace the dispersion in the vicinity of Ep. The 
energy gap was determined from the difference between 
Ep and the maximal energy that a back-bending disper- 
sion reaches [Fig.[4](g)]. Fig.[4^i) shows the results of such 
an analysis applied to the data taken at 54 K. Within the 
error bars, the obtained gap function is consistent with 
that derived from the symmetrization method. 

To explain our results, one natural line of thinking is to 
assume that both the superconducting gap and the pseu- 
dogap have d x 2_ y 2 form, and the combination of disor- 



3 



(a) 

137 K 
88 K 



(b) 




■x = 0.08, 10 K 

x = 0.08, 54 K 

x = 0.03, 10 K 



-0.12 -0.08 -0.04 0.00 0.04). 12 -0.08 -0.04 0.00 0.04 
energy (eV) energy (eV) 




-0.12 -0.08 -0.04 0.00 0.040.12 -0.08 -0.04 0.00 0.04 
energy (eV) energy (eV) 





0.28 0.30 0.32 0.34 0.3 
k (jt/a) 



-0.20 -0.15 -0.10 -0.05 0.00 
energy (eV) 



FIG. 2: (Color online) EDCs at Jzf on the zone diagonal 
and band dispersions along the zone diagonal for LSCO. The 
EDCs were obtained by deconvolving the raw ARPES data 
to remove the broadening due to the finite instrumental res- 
olution 4 ]. (a) EDCs at /cf as a function of temperature for 
x — 0.08. Curves are offset vertically for clarity, (b) EDCs 
at k F for x = 0.03 at 10 K, and for x = 0.08 at 10 K and 54 
K. Curves are offset horizontally to align the peak position to 
that of the EDC for x = 0.08 at 10 K. (c) Dispersions derived 
from the peak positions of the momentum distribution curves 
for x = 0.08 at 10 K (< T c ) and 54 K (> T c ). The kink 
is indicated by the black arrow, (d)-(e) EDCs at 10 K as a 
function of doping (x). Curves in (d) are offset vertically for 
clarity and in (e) are offset horizontally to align the peak po- 
sition to that of the EDC for x = 0.08 at 10 K. (f ) The same 
EDCs as in (d), but the spectra are normalized to intensities 
at high binding energy. 



der with long-range Coulomb interactions depresses the 
density of single-particle excitations at Ep, forming a so- 
called Coulomb gap [12j [13] . A counter argument to this 
picture is that the scattering due to the disorder should 
also reduce the lifetime of quasiparticles, which would 
be manifested as severe broadening in the ARPES line- 
shapes. Such broadening is not seen in the data, where 
instead the widths of the EDCs are the same in the su- 
perconducting state for samples with different dopings. 
Moreover, these widths are smaller than those at T > T c 
or in the non-superconducting sample at low temperature 
(Fig. [2j. One could instead consider a scenario in which 
an ordered phase (e.g., a spin-density wave [HJ [15]) or 
its flutuations above T c opens an energy gap along the 
entire FS in lightly doped LSCO (x < 0.105). In such a 
picture, this phase is different from the superconducting 
instability, and below T c it coexists with the supercon- 
ducting phase, which has a d x i_ y i order parameter. The 
end form or the combination of these two order parame- 
ters would produce an energy gap on the entire FS. Re- 
cent Monte Carlo simulations have shown that a strong 
fluctuating competing order could open a diagonal gap 
in the superconducting phase, and a d x i_ y i gap function 
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FIG. 3: (Color online) ARPES spectra for LSCO (x = 0.08). 
(a) - (c) Symmetrized EDCs as a function of temperature 
at k F (a), (b) and (c) shown in (d). (d) The first quad- 
rant of the BZ. The red curve is the FS. The closed circles 
indicate the kF at which the symmetrized EDCs are shown 
in (a)-(c). (e) Experimentally determined energy gap as a 
function of the FS angle <j) indicated in (d) at 10 K (circles) 
and 89 K (diamonds). The dash line is the d x i_ y i form in 
38.52 meV. The dotted line is the 
with A2 = 20 meV. The solid line is 
I Ad 2 _ 2 + 77" (f) Energy gap as a function of <j). The 
symbols are the experimentally determined energy gaps at 
various temperatures. The solid lines are |A<j 2 _ 2 + ^d xy \ 
with A° d = 38.52 meV and A^ = 20 meV (red), 16.6 
meV (green), 13.5 meV (blue), and meV (yellow-green). 
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is restored with reduced amplitudes of the fluctuating 
competing order [16 . This scenario would explain why 
the diagonal gap is open in LSCO x = 0.08 and 0.03 at 
low temperatures and is closed for x > 0.1, since it has 
been suggested that disordered magnetism is present at 
low doping and disappears at optimal doping [17j ITS] . 
The open question is how to account for the similar mo- 
mentum dependence of the gap above and below T c in 
LSCO (x — 0.08), as well as the fact that the ampli- 
tude of the gap at the anti-node is relatively insensitive 
to temperature. An alternative explanation is that in 
LSCO there is a critical doping point below which the 
superconducting gap function changes from d x i_ y i to an- 
other form. The pseudogap at T > T c is a precursor 
to superconductivity [19-22 , which has similar momen- 
tum dependence and amplitude as the superconducting 
gap. To study the symmetry of the gap function in the 
superconducting state of the highly underdoped LSCO 
(x = 0.08) sample, we have compared the obtained gap 
function at 10 K [Fig. [3]^e)-(f)] with those of all allowed 
spin-singlet even-parity pair states in cuprates [3] [23]. 
The only order parameter basis function that is nodeless 
in its pure form is the s- wave state. However, although 
this order parameter is nodeless, its momentum depen- 
dence is not compatible with the observed momentum 
dependence of the energy gap. A combination of two 
real subcomponents of different superconducting order 
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FIG. 4: (Color online) ARPES spectra of 
LSCO (x = 0.08) measured at 54 K (T c = 
20 K). (a)-(f) Intensities along 6 cuts in- 
dicated in (h). The spectra were obtained 
by the DFD method, (g) Dispersions ob- 
tained by tracing the EDC peaks of (a)- 
(f) in the vicinity of Uf. For clarity the 
dispersions along cuts 2-6 were offset hor- 
izontally, (h) The FS (red line) obtained 
from TB fits to the experimental data. 
Lines 1-6 indicate the cuts along which 
the data in (a)-(f) were taken, (i) Energy 
gap as a function of FS angle (</>). The 
closed circles are the gap extracted from 
the difference between Ef and the maxi- 
mal energy that the back-bending disper- 
sions reach. The triangles are the gap ob- 
tained from symmetrized EDCs [Fig. |3jf)]. 
The solid lines are the | Ad x2 _ 2 | (blue) and 
l A ^2_„ 2 + iA ^J (green) as in Fig. [3^). 



parameters, such as a mixed s + d x 2_ y 2 gap [24 , will 
only shift the node(s) along the FS, which cannot pro- 
duce a nodeless gap structure. By considering all the 
possible nodeless mixed gap function we find that the 
combination of d x 2_ y 2 + id xy gaps is the only form that 
can reproduce the momentum dependence of the mea- 
sured gap function quantitatively. In Fig. [3je) we plot 
| Ad 2 _ 2 (k)\ and |A^ (fe)|, as well as the modulus of 
their sum A(fe) = \Ad 2 _ 2 (fe) + iAd xy (fe)|, as a function 
of FS angle (</>), where 

^d xy (fe) = A2 av [sm(k x a) sin(V)] • (2) 

We take A[] 2 = 1.07A antinode and A^ = A diag , 
where A antino d e and Adiag are the measured gap sizes 
at the zone boundary and the zone diagonal, respec- 
tively. The pre-factor of 1.07 comes from 2/ [cos (kp x a ) ~ 
cos(kFyd)] where (kFxi^Fy) is the intersection of the FS 
and the zone boundary. Remarkably, excellent agreement 
is found between the experimental data and the mixed 
d x 2_ y 2 +id xy gap function along the entire FS for all the 
temperatures. To see this, we simply hold A 9 con- 

— y 2 

stant and use the measured gap sizes on the zone diagonal 
(Adiag) for the amplitude of the d xy contribution to the 
mixed d x 2_ y 2 + id xy gap at different temperatures. The 
results are shown in comparison to the data in Fig. [3]^f). 
Interestingly, the mixed d x 2_ y 2 + id xy pair state implies 
that the time-reversal symmetry is broken in the system 
[3], which has been embodied in early theories of quan- 
tum phase transitions [25j [26] . 

To summarize, our main experimental findings are: 1) 
for highly underdoped LSCO (x = 0.08), in the supercon- 
ducting state the electronic excitations are gaped along 
the entire underlying FS; 2) the diagonal gap persists 
above T c while increasing the temperature and/or reduc- 
ing the hole-concentration; 3) whereas the gap size on 



zone boundary remains constant up to 150 K, the diago- 
nal gap is temperature-dependent; 4) when the diagonal 
gap decreases to zero, it leaves behind a pure d x 2_ y 2 en- 
ergy gap, at which point a Fermi arc emerges and its 
length increases with temperature. Our observations in 
highly underdoped LSCO could be explained either by 
a strong fluctuating competing order which is different 
from the superconducting order parameter [16 , or by a 
mixed d x 2_ y 2 + id xy gap function occurring when the 
doping is below a quantum critical point [25 j. The en- 
tirely gapped Fermi surface in highly underdoped LSCO 
is a profound departure from the d x 2_ y 2 form observed 
at moderate to high doping. Further investigations of the 
origins of this nodeless gap function and its evolution to 
pure d x 2_ y 2 form promise to shed light on how super- 
conductivity emerges from the Mott insulating state in 
high-T c cuprates. 
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